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Various spin effects are expected to become observable in light-matter interaction at relativistic intensities. 
Relativistic quantum mechanics equipped with a suitable relativistic spin operator forms the theoretical foundation 
for describing these effects. Various proposals for relativistic spin operators have been offered by different authors, 
which are presented in a unified way. As a result of the operators’ mathematical properties only the Foldy- 
Wouthuysen operator and the Pryce operator qualify as possible proper relativistic spin operators. The ground 
states of highly charged hydrogen-like ions can be utilized to identify a legitimate relativistic spin operator 
experimentally. Subsequently, the Foldy-Wouthuysen spin operator is employed to study electron-spin precession 
in high-intensity standing light waves with elliptical polarization. For a correct theoretical description of the 
predicted electron-spin precession relativistic effects due to the spin angular momentum of the electromagnetic 
wave has to be taken into account even in the limit of low intensities. 


1. Introduction 

Employing novel light sources such as the ELI-Ultra High Eield Eacility, for example, that envisage to provide field intensities 
in excess of 10^*' W/cm^ and field frequencies in the x-ray domain HHU light-matter interaction in the relativistic regime may 
be probed experimentally. Relativistic quantum mechanics predicts various new phenomena to occur in this regime El El, for 
example, multiphoton scattering, radiation reaction effects, vacuum-polarization effects or even pair creation EHH. Furthermore, 
electrons in strong electromagnetic fields can exhibit distinct spin effects MM- 

Relativistic quantum mechanics has to be employed to study spin phenomena in strong electromagnetic fields. According to the 
formalism of quantum mechanics, each measurable quantity—as the spin, for example—is represented by a Hermitian operator. 
However, there is no universally accepted operator to describe an electron’s spin degree of freedom within the framework of 
relativistic quantum mechanics. Thus, we first investigate the properties of different proposals for a relativistic spin operator and 
show that most candidates are lacking essential features of proper angular momentum operators IfTTlfTSlI . Only the so-called 
Foldy-Wouthuysen and the Pryce operators qualify as proper relativistic spin operators. The various spin operators predict 
different expectation values when electrons interact with electromagnetic potentials. In this way, one may distinguish between the 
proposed relativistic spin operators by experimental means. In particular, eigenstates of highly charged hydrogen-like ions may be 
utilized to identify a legitimate relativistic spin operator experimentally. 

A further relativistic spin phenomenon, which we study in more detail, is the coupling of the spin angular momentum of light 
beams with elliptical polarization to the spin degree of freedom of free electrons IT9ll20l . This coupling, which is of similar origin 
as the well-known spin-orbit coupling, and the magnetic field lead to electron-spin precession. The spin-precession frequency 
is proportional to the product of the laser-field’s intensity and its spin density. To derive the correct spin-precession frequency 
relativistic corrections to the nonrelativistic Pauli equation, which account for the light’s spin density, have to be taken into 
account. The quantum mechanical interactions of the electron’s spin with the laser’s rotating magnetic field, which may be 
characterized by the nonrelativistic Pauli equation, and the electron spin’s interaction with the laser field’s spin density, which 
results via the relativistic corrections, counteract each other. As a result, a net electron-spin rotation remains with a precession 
frequency that is much smaller than the frequency predicted by a nonrelativistic theory. These relativistic effects are maintained 
even if the involved electromagnetic field strengths are nonrelativistic. 


2. Relativistic spin operators 

A Lorentz invariant quantum mechanical description of the motion of an electron in electromagnetic fields is provided by the 
time-dependent Dirac equation. For a particle of rest mass mo and charge q it is given by (units are used in this section for which 
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n^i) 


= H'V{r, f) = (cor ■ (p - qA(r, t)) + q(f>ir, t) + nioc^p) 'P(r, f) , (1) 

with the electromagnetic potentials (p{r, t) and A{r, f), the speed of light c, the canonical momentum operator p - -iV, and the 
matrices a - {ai,a 2 , and (3. These 4x4 matrices obey the algebra 

= 1 , atOk + a^ai - 25 ,- ^, afi + [3ai = 0 . ( 2 ) 


To specify our notation and abbreviations, we employ the Dirac representation, where the matrices a,- and p are defined as 


0 (Ti 
o-i 0 


for / = 1 , 2 ,3 , p 


I 2 0 

0 -I 2 


in terms of the three 2x2 Pauli matrices tr = (crj, 0 - 2 , 0 - 3 )'. In the Dirac representation, the Pauli matrices are given by 


CTl 


0 1 

1 or 


0-2 


0 -i 
i 0 ’ 


0-3 


1 0 
0 -1/ ■ 


(3) 


(4) 


The symbol I 2 denotes the 2x2 identity matrix. The free-particle Dirac Hamiltonian with A(r, t) — 0 and 0(r, t) — 0 will be 
denoted by Hq. The free-particle Dirac Hamiltonian features positive as well as negative energy-eigenvalues. We will also use the 
operator po to denote 


PO = (5) 

which has the same eigenstates as Hq but all its eigenvalues being positive. Furthermore, the matrices 

^'"( 0 ' fori =1,2,3 ( 6 ) 

will be employed, which will be commonly combined into the three-component operator S = (Si, S 2 , 2 ^ 3 )^. 

A relativistic spin operator may be introduced by splitting the undisputed total angular momentum operator J, which is in the 
Dirac representation 


/ = rx(-iV)+is, (7) 

into an external part L and an internal part S, viz., J - L + S. These parts are commonly referred to as the orbital angular 
momentum and the spin. Because the orbital angular momentum operator L is related to the position operator f and the momentum 
operator p — -iV via L - fx p, different definitions of the spin operator S imply different relativistic position operators f. Thus, 
the question for the right splitting of the total angular momentum into an orbital part and a spin part is closely related to the quest 
for the right relativistic position operator 

There are two complementary approaches to determine a suitable relativistic spin operator from a set of possible candidates. A 
common approach is to judge a spin operator candidate by its mathematical properties. This means to analyze its symmetries, 
its behavior under various transformations, its algebraic features, and so on. Historically, most relativistic spin operators have 
been proposed on purely mathematical arguments. Considering that relativistic spin effects are expected to be detectable at high- 
intensity laser facilities one may also compare theoretical predictions based on various relativistic spin operators to experimental 
measurements. In this way, candidates for relativistic spin operators that are not compatible with measurements can be ruled out 
on the basis of experimental results, rather than solely by mathematical reasoning. Here, we will adopt both approaches. 

Mathematically, we demand from a proper relativistic spin operator S - (5 1 ,5 2 ,5 3 )^ and its three components the following 
features; 

1. Each component of a spin operator should be a Hermitian operator. 

2. The physical quantity that is represented by the operator S should not depend on the orientation of the chosen coordinate 
system. Thus, a spin operator must transform under rotations like a vector, which is ensured by fulfilling Il50l 

[Ji,S p^iSkj^kSk (8) 


with Sij,k denoting the Levi-Civita symbol. 
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3. It is also required to commute with the free Dirac Hamiltonian, i. e., [Ho,S] = 0. This property ensures that the relativistic spin 
operator is a constant of motion if forces are absent, such that spurious Zitterbewegung of the spin is prevented. 

4. A spin operator must feature the two eigenvalues +1/2 and it has to obey the angular momentum algebra 

[5i,5j] = ieij,A:5*. (9) 

These two requirements are commonly regarded as the fundamental properties of angular momentum operators of spin one-half 
particles ISOl . 

Table[T]gives an overview over several spin operators, which have been proposed in the literature, and presents their mathematical 
key features. A detailed description of these operators and their relations to each other is given elsewhere ifTSlI . Note that some 
of these spin operators have been discovered by different authors in different contexts yielding different but mathematically 
equivalent forms. For example, the so-called Newton-Wigner spin operator ED, which may be written as ED 


A _ Po ^ 
•JFW - ^ 


pi. 


IniQC 2moc{moc + po) 


p- pxa 


iHo 


Inioc^po ’ 


( 10 ) 


is just another way to express the Foldy-Wouthuysen spin operator given in Tab. 

The Foldy-Wouthuysen, the Chakrabarti, the Pryce, and the Fradkin-Good spin operators are equivalent in the positive-energy 
subspace of free-particle states. This can be shown easily in the Foldy-Wouthuysen representation. The transition from the 
standard representation, where the Dirac Hamiltonian has the form Q and the various spin operators have the forms given in 
Tab.[T] is mediated via the nonlocal unitary transform 


Tfw — 


Po + nioc - j3a ■ p 
^J2po(po + moc) 


In the Foldy-Wouthuysen representation the free-particle Dirac Hamiltonian becomes diagonal, 

^0 “ ^fw^oTfw = cPpQ 


and the Foldy-Wouthuysen spin operator has the form 


■^FW “ ^p^SfwTfw - 


jtr 0 

I 0 lo-. 


( 11 ) 


( 12 ) 


(13) 


Since the momentum operator is invariant under the Foldy-Wouthuysen transformation fpw, the simultaneous eigenstates of the 
free-particle Dirac Hamiltonian, the momentum operator, and the z component of the Foldy-Wouthuysen spin operator are in the 
Foldy-Wouthuysen representation given by 



Ti 


(Oi 


(Oi 


(Of 

*FW,+,p,T - 

0 

0 

2’^ ^ , *FW,+,p,i - 

1 

0 

Sfw,-,p,t - 

0 

1 

'', *FW,-,p,i - 

0 

0 


loj 


loj 


loj 


UJ 


where the indices indicate the sign of the energy eigenvalue, the momentum eigenvalues and the spin eigenvalue (t for 1 /2, 
I for -1/2). In the Foldy-Wouthuysen representation the Chakrabarti, the Pryce, and the Fradkin-Good spin operators are given 
by the rather simple expressions 


■^Ch “ ^fW’^ChTfw - 


Spi- - fp^Spr^FW - 


/ 1 • <TXp 

k(T 1 —- 

2 niQC 

,0 io-j’ 
0 

0 


- I - - ^ 

S f _ 'T'—l C _ 4, 

FG “ ^fW‘^FG^FW “ 


0 


\(Tj 


(15) 

(16) 
(17) 


Because the operators ([T^, ( [TSl l, ( [Tbl l, and ( [TtI i have the same upper left 2x2 matrix, all these operators act in the same way on 
the positive-energy free particle states given in For quantum states that are superpositions of free-particle states with positive 
and negative energy, however, these spin operators are not equivalent. 
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On the basis of the four criteria given above, one may argue that only the Foldy-Wouthuysen spin operator and the Pryce spin 
operator qualify as proper relativistic spin operators because only these two fulfill all four criteria. However, the question of 
which of the proposed relativistic spin operators in Tab. [^provides the correct mathematical description of spin can be answered 
definitely only by comparing theoretical predictions with experimental results. Because if an interaction with some external 
fields is introduced a superposition of positive-energy free-particle states evolves such that negative-energy free-particle states 
become populated. The proposed spin operators are not equivalent if they are applied to positive-energy states of Hamiltonians 
with nonvanishing electromagnetic fields. Therefore, it becomes possible to distinguish between the various spin operators by 
determining their expectation values for electrons interacting with electromagnetic fields. 

For this purpose it is desirable to employ a physical system that shows strong relativistic effects and is as simple as possible. 
Such a setup is provided by the bound eigenstates of highly charged hydrogen-like ions, i. e., atomic systems with an atomic core 
of Z protons and a single electronic charge. These ions can be produced at storage rings Il52ll or by utilizing electron beam ion 
traps Il5^l54l up to Z = 92 (hydrogen-like uranium). The degenerate bound eigenstates of the corresponding Coulomb-Dirac 
Hamiltonian (in atomic units) 


He ^Hq-- 


( 18 ) 


are commonly expressed as simultaneous eigenstates of He, J^, Jj,, and the so-called spin-orbit operator K - 

f3{t, ■ [r X (-iV) + 1)]) fulfilling the eigenequations Il55ll5^ 


Hc4^n,Kj,m — ^(n, n,K,j,m 5 

J — jij "t" i 

The eigenenergies are given with cTei denoting the fine structure constant by 


Sin, j) = moc^ 


1 -H 




n = 1,2,... , 

|/c| = 1 , 2 ,..., n, /<■ —n , 

; = kl - 5, 

'« = - j, ij - 1), ■ • ■,; ■ 

- 1/2 


1/2+ VO'-1/2)2-a2z2 

The degenerate hydrogenic ground state is with y = Vl “ radial function 


^-moZr 


Hr) 


{2mi)Zry r ’ 


(19a) 

(19b) 

(19c) 

(19d) 


( 20 ) 


( 21 ) 


and the normalizing factor 


given by the two wave functions El 


TV = (2moZ)2^2 


1 -I- 7 


(2F(1 -H2y) 


( 22 ) 


ki 1 i i{r,0,4)) = NHr) 


ki,i i= TVi/r(r) 


YoM4>) ^ 

0 

i^V^Fi.o(k<^) 

0 

YoMH 

i^ 4 F,_,(e, 0 ) 


(23a) 


(23b) 


The spin expectation values of the z component of the in Tab. [T] defined spin operators are presented in Fig.[2a). For small 
atomic numbers (Z % 20), all spin operators yield about 1 /2; for larger Z when relativistic effects set in, however, expectation 
values differ significantly from each other. While for Pauli, Fouldy-Wouthuysen, Czachor, Chakrabarti, and Fradkin-Good spin 
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Figure 1 : (a) Spin expectation values of various relativistic spin operators for the hydrogenic ground state \23n) as a function of the atomic 
number Z, adopted from New J. Phys. 16(4), 043012 (2014). For the ground state ( |23b| ) we find the same spin expectation values but with 
opposite sign (not displayed in the plot), (b) Expectation values of the position variance with respect to the z coordinate for the hydrogenic 
ground state p3a| > as a function of the atomic number Z. The wave function’s variance depends on the definition of the relativistic position 
operator, which is induced by the definition of the relativistic spin operator. 


operators the spin expectation value is reduced, the expectation value of the Frenkel spin operator exceeds 1/2. Only for the 
Pryce operator we find that the spin expectation values is 1 /2 for all values of Z, which can be also shown analytically IfTFI . A 
comparison of the results of a spin-measurement experiment for groundstate electrons in hydrogen-like highly-charged ions to the 
numerical results Fig.[T]a) would allow to find a suitable relativistic spin operator or at least to rule out some candidates ifTTl . 

As outlined above, each definition of a relativistic spin operator S induces also a relativistic position operator r. Because only 
the Fouldy-Wouthuysen and the Pryce operators fulfill all our mathematical criteria for a proper spin operator we focus on these 
two in the following. For determining the corresponding position operator the defining relation 

7 = r X (-iV) + 1112 ^rx (-iV) -h S (24) 


is rather unwieldy. In the case of the Fouldy-Wouthuysen and the Pryce operators it is more convenient to utilize the fact that 
these are related to the Pauli spin operator via a unitary transform. The position operator of the Pauli spin operator is just 


rp^r. 

Therefore, the position operators of the Fouldy-Wouthuysen and the Pryce operators are given by 

= TFwfTp-^ 

and 

rpr = fprrfpj' 

with the transformations Ffw and Tpr defined in ( [TT] i and by 

(h 0 ^ 

fp.= ^ 

m 


(25) 


(26) 


(27) 


(28) 


Determining an explicit expression for the transformation (|26ll yields the famous Fouldy-Wouthuysen mean position operator 


rpw - r + ih. 


i£ X p 


Pia ■ p)p pa ^ 


[2po(po + niQc) 2pl(po + moc) 2po 

which is equivalent to the Newton-Wigner position operator ll2Tll . while ^n\ yields 

fO 0 

rp^ - r -fi 


(29) 


0 


ax p 

IpP 


(30) 
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For numerical calculations of expectation values, however, it is more convenient to utilize the forms ( |26] l and ( |27] i. 

Because of symmetry reasons the expectation values of the position operators rp, rpw, and rp, are zero if these position 
operators are applied to the hydrogenic ground states ( |23a[ l and ( |23b[ ). The second moment, however, does not vanish. The 
variance of a nonrelativistic hydrogenic ground state scales with 1 /Z"^. Due to relativistic effects the ground states’ wavepackets 
shrink even faster. How fast, depends on the definition of the position operator, as shown Fig.fTlb). In principle, one can measure 
the variance of the ground state of a highly charged ion and compare it to the predictions in Fig^^Jb) and in this way determine the 
correct relativistic position operator and indirectly also the correct relativistic spin operator. 


3. Electron-spin precession in elliptically polarized light 

In the previous section we argued that the spin of free electrons should be modeled by the Fouldy-Wouthuysen spin operator. In 
the following we will utilize this operator to study a relativistic spin dynamics, which originates from a coupling of the electron’s 
spin to the spin of an electromagnetic wave with elliptical polarization. 

The electric and magnetic field components of two elliptically polarized laser fields propagating into the positive or negative 
direction of the x axis are given by 


y N 1 27r(x + ct) 

Ei, 2 {r,t) = £ cos-;-Cv + cos 


= — l + cos 


A 

27t{x + ct) 


2’k{x + ct) 


±rj\e 


2n{x + ct) 

+ 7/1 Cy + COS -C; 


(31a) 

(31b) 


Here, the position vector r - {x,y,z)^, the time f, and e^, Cy, and denoting unit vectors in the direction of the coordinate axes 
are used. The parameter rj e (-tt, tt] determines the degree of the light beams’ ellipticity with t] - 0 and rj = tt corresponding to 
linear polarization and rj - +7r/2 to circular polarization. The two electromagnetic waves ( [3T] i feature the same wavelength A, the 
same electric field amplitude E, and the same intensity 

I = eocE^ , (32) 


but have opposite helicity. Introducing the wave number k — 2nlA and the lasers fields’ angular frequency cu = kc, the Coulomb 
gauge vector potentials Ai 2 (r, t) of the elliptically polarized fields pT|) are 


E / X 

Ai 2 (r, t) - -(+ sin(A:x + 0 Jt)ey + sin(A:x + cot ± r])eA . 

Each of the electromagnetic fields specified by carries the photonic spin density 

_ ^ soE^A^inri 

X ^ 1,2 - -;;-• 

Inc 


(33) 


(34) 


As one can show via the Volkov solution of the Dirac equation |59ll a single plane wave as given in ( |3T] i cannot change the spin 
orientation of an electron. Therefore, we consider a standing wave, which is formed by superimposing the two counterpropagating 
waves given in (|3T]i. The magnetic vector potential of the combined laser fields is given by 


A(r,f) = - 


2w{t)E 


i kx ^sin (jjt By + sin(a;f - rj) . 


(35) 


Here the window function 


ifO < f < AT, 

if AT < f < r - AT, (36) 

\fT - AT <t <T, 

was introduced to allow for a smooth turn-on and turn-off of the laser field. The parameters T and AT denote the total interaction 
time and the turn-on and turn-off intervals. For circularly polarized plane waves ir] - 7r/2), the electric and the magnetic 
components of the standing wave are parallel to each other and rotate around the propagation direction. The maxima of the 
electric and the magnetic field components are shifted against eachother by AjA. 

Solving the time-dependent Dirac equation till time t — T with a common eigenstate of the free Dirac Hamiltonian, the 
momentum operator, and the z-component of the Foldy-Wouthuysen spin operator with zero momentum and positive spin as 


wit) = 


Sin*' 

1 

sin 


2^T 


2 n{T~t) 
2 Ar 











8 


10 *' 


10 ^ 


I (W/cm^) 


10 "* 


10 *' 


10 “ 


G jq13 


10 “ 


10 ‘ 


O Dirac eq. 

+ relativistic Pauli eq. 
□ Pauli eq. 


10 “ 
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Figure 2: Angular frequency of the spin precession as a function of the laser’s electric field strength E and its intensity / for electromagnetic 
fields with a wavelength A = 0.159nm. Depending on the applied theory (the Dirac equation Q. the relativistic Pauli equation l |38[ l, or the 
nonrelativistic Pauli equation) the spin of an electron in two counterpropagating circularly polarized light waves scales with the second (Pauli 
equation) or the fourth power (Dirac equation, relativistic Pauli equation) of E. Numerical data adopted from New J. Phys. 16(4), 043012 (2014). 


initial condition shows that the electron’s spin precesses around the propagation axis of the electromagnetic fields lfT9ll20l . The 
role of the photonic spin density for the electronic spin precession becomes evident by considering the weakly relativistic limit of 
the Dirac equation Q. In this limit, this equation reduces via a Foldy-Wouthuysen transformation lIMlISTIIhOl to 


iH'¥(r, t) = 


(-iW - qA{r, t)y qh 


2mo 


- — o- ■ B(r, t) + q(f>{r, t) - 
IniQ 


(-im - qA(r, f)r 


^ ^ (c^B(r, tf - E(r, tf) 




o- ■ {E{r, t) X (-iliV - qA{r, t))) - ■ E(r, t) + 






qti 


{o- ■ B(r, t), i-im - qA{r, t)f] T'(r, t) (37) 


for the now two-component wave function 'I'(r, t) with the vector potential A{r, t) given by ( |35] l and the electromagnetic fields 
B(r, f) = V X A(r, f) and E{r, f) — -A{r, t). In leading order Eq. ( [JT) ! features four terms that may cause spin dynamics. The 
so-called Zeeman term ~ cr ■ B(r, t) and its the lowest-order relativistic correction given by the anticommutator expression 
~ {(T ■ B{r, t), (-i^V - qA(r, f))^) mediate the coupling of the electron’s spin to the magnetic field. The term ~ cr ■ E(r, t) x i/iV 
leads to the so-called spin-orbit interaction, i. e., the coupling between the electron’s spin and its orbital angular momentum. 
The term ~ cr ■ {E{r, t) x A(r, f)) may be interpreted as a coupling of the spin density of the external electromagnetic wave to 
the electron’s spin. Considering the relativistic correction due to the electromagnetic wave’s spin density as the only relativistic 
correction to the nonrelativistic Pauli equation the relativistic electron motion in the vector potential ( |35] ) may be described by the 
relativistic Pauli equation ifT^ l20l 


imiir, t) = (^(-iW - qA{r, t)f - ■ B{r, t) + ^ cr ■ {E{r, t) x A(r, f))] T'(r, t). (38) 

\2mo 2mo / 

Numerical calculations indicate that the relativistic Pauli equation ( [38] l is sufficient to reproduce the spin dynamics of the fully 
relativistic Dirac equation Q. The nonrelativistic Pauli equation, however, yields a completely different spin dynamics as shown 
in Fig. 1^ Depending on if the Dirac equation Q, the relativistic Pauli equation ( [38l l, or the nonrelativistic Pauli equation is 
applied, the spin of an electron in two counterpropagating circularly polarized light waves scales with the second (Pauli equation) 
or the fourth power (Dirac equation, relativistic Pauli equation) of the electric field’s amplitude E. Further analytical calculations 
based on the Dirac equation and time-dependent perturbation theory ll20l show that the electron’s spin precesses with an angular 
frequency that is proportional to the photonic spin density Qo- - eoE^Ajinc), the laser field’s intensity I given in ( [3^ , and the 
fourth power of the wavelength: 


Q = Po-// 




( 39 ) 


The proportionality factor a^j/(27r^mgC^) is independent of the standing wave’s electromagnetic field. 











9 




Figure 3: An ensemble of classical spins under the effect of a magnetic field and a photonic spin density of a standing wave formed by two 
counterpropagating circularly polarized waves. Part a): Initially spins are aligned. Part b): The rotating magnetic field causes a position- 
dependent rotation of the spin vectors. Part c): Also the photonic spin density causes a spin rotation but into the opposite direction. Averaged 
over a wavelength both effects cancel eachother. 


The role of the relativistic correction due to the photonic spin density is pivotal, its influence on the scaling of the spin 
precession frequency does not become small in the limit of weak fields, which is usually related to a nonrelativistic limit. The 
nonexistence of a nonrelativistic limit results because the Zeeman term is due to the fast oscillation of the magnetic held effectively 
as strong as the relativistic correction due to the photonic spin density. On the basis of a classical argument lfT9ll one can show that 
the effects of the Zeeman term and the correction due to the photonic spin density on the motion of the electron spin counteract 
eachother. The initial electron’s quantum state, which is a momentum eigenstate, is delocalized over several laser wavelengths. 
Thus, we mimic the quantum wavepacket by an ensemble of classical particles with spin angular momentum. These particles are 
placed along the x axis with initially aligned spin direction, see also Fig.j^a). The dynamics of a classical electron spin s at fixed 
position r in the magnetic held B{r, t) - Bi(r, t) + B 2 {r, t) is governed by the classical equation of motion 


s{t) = —s(t) X B(r, t) = 
mo 


2qE cos kx 

- s{t) X (cos cot Cy + sm cotCy). 

mo 


(40) 


In a static and homogeneous magnetic held of strength E/c, the electron’s spin would precess around an axis parallel to the 
magnetic field’s direction with the angular frequency Ql = qEf{moc). For a rotating magnetic held and parameters such that 
Ql « CO, however, the magnetic held rotates so fast that the spin precesses around the rotation axis of the magnetic held. For the 
magnetic held of a standing wave formed by two counterpropagating circularly polarized waves given by the vector potential ( |T5] ) 
a position-dependent angular frequency 2Qp sin^ kx results with Op = {qE)^Aj see Fig.|^). Similarly, the relativistic 
correction due to the photonic spin density in ( [38| l to the Pauli equation leads to the classical equation 

cP' {qEpA cos^ kx 

X {E{r, t) X Air, t)) = ’ -s(f) x e,. (41) 

Zm^c^ TTm^c^ 

This yields also a position-dependent angular frequency, but now -2Qp cos^ kx, see Fig.j^c). As a consequence, a classical spin 
under the effect of both spin terms, the Zeeman term and the relativistic term due to the photonic spin density, rotates in a short 
time interval At around an angle of about 2Qp(sin^ kx - cos^ kx)At. Averaged over a laser wavelength this rotation angle vanishes 
and the effects of the two spin terms chancel each other in our classical model. Thus, the classical model explains how the effect 
of the laser fields’ spin density leads to a breakdown of the quadratic scaling of the spin-precession angular frequency in E that 
results if only the magnetic held is taken into account. The model, however, is not able to reproduce the quartic scaling in E that 
results from the fully relativistic quantum mechanical Dirac equation. The quartic scaling results as a genuine quanmm effect 
from the fast temporal oscillations combined with the spatial modulation of the electromagnetic fields. In fact, further numerical 
calculations show that the quantum mechanical wavepacket accumulates in regions of the standing light wave with high magnetic 
fields 1 ^ . 
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